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1. Introduction 


The analysis of neutrino oscillation experiments requires having under control all the nuclear 
effects, which are inherent to any v-nucleus scattering event. These effects play a major role in 
modifying the free v-nucleon cross section, and a good understanding of them is clearly mandatory 
in order to reduce the systematic uncertainty of the models employed in the experimental analysis. 
There is strong theoretical evidence [Q, about a significant contribution from multi-nucleon 
knockout to the inclusive Charged Current (CC) cross section around and above the quasielastic 
(QE) peak region. 

There are, at least, three different microscopic models which are based in the relativistic 

Fermi Gas. But they differ from each other in several assumptions and different nuclear ingredients 
for the interaction. Therefore, it is really difficult to disentangle the origin of any discrepancy in 
the final resulfs. 

Ofher models 


101 are based on fhe available phase space Jusf assuming a consfanf fransifion 


mafrix elemenf and filling if lo fhe experimenlal cross section. 

Finally, our goal is lo reduce fhe compulalional time needed in Ihis kind of calculalions, or al leas! 
lo eslablish whal assumptions previously made by olher aulhors are really good enough in order 
lo estimate accurately and fasl Ihe conlribulion of Ihese mulli-nucleon processes lo Ihe inclusive 
channel. 


2. 2p-2h phase space in the Relativistic Fermi Gas model 


The hadron tensor for Ihe 2p-2h channel in a fully relativistic framework is given by: 


^''2p2h 


I d^p\ d^p'^ d^h2 (P'l ,v’2MM)5{E[+E'^-Ei-E2-co) 

+P 2 -hi -h 2 -q) (2.1) 


where is Ihe nucleon mass, V is Ihe volume of Ihe system and we have defined Ihe producl of 
step functions, 

&{p\,P2,hi,h2) = G{p[-kf)e{p2-kF)6{kF -hi)e{kf -h2) (2.2) 


which encodes Ihe nuclear model. 

Finally, Ihe function r^''(pj,p 2 ,hi,h 2 ) is Ihe elemenlary “hadron” tensor for Ihe fransifion of a 
nucleon pair wilh given initial (hi,h 2 ) and final (pj,p 2 ) momenfa, summed up over spin and 
isospin, given schemalically in terms of fhe anlisymmelrized fwo-body currenfs by: 


Kt''' 


(p'i,P2,hi,h2) 


^£/(1',2';1,2);/(1',2';1,2),4 

' rr 'T 


(2.3) 


The above multidimensional inlegral (^ri]) can be done eilher numerically or ifs dimensions can be 
furlher reduced under some approximations SOl-IIl- We do no! know exaclly fhe origin of fhe 
discrepancies befween fhe available models. These could be due fo differenl fwo-body cuirenfs, fo 
local Fermi Gas model rafher fhan global one, or ofher effecfs. If is difficulf fo make any concluding 
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assessment right now, but all those models should agree at the level of phase space function F(ft), ^), 
obtained assuming a constant 


F{( 0 ,q)= J d^hid^h 2 d^ p\ ^-^^j^@{p\,P2,hi,h2)5{E[ + E 2 


-E 1 -E 2 -C 0 ) (2.4) 


where = 1 and P 2 = hi + h 2 + q — integrating out the delta function of momentum conser¬ 
vation. 

The remaining delta function enables analytical integration over the modulus of p^: 


E{(0,q)=27i I d^hid^h 2 de[smd[^^ ^ 1 

J a=±F2Pl ■^iP2'Pll 


2 12 
mffPi 


(2.5) 


/ '(“) 
Pi=Pi 


and the sum inside the integral sign runs over the two solutions of the energy conservation 
delta function (see appendix C on Ref. [[T]]). 


2.1 Frozen nucleon approximation 

The frozen nucleon approximation is just a particular case of the mean-value theorem in several 
variables 



f{c){b — a) with c£[a,b\ 


f /(r)ri"r = /(c) f ri”r = /(c)r 
Jt Jr 


with c G y 


( 2 . 6 ) 

(2.7) 


In our case we are going to skip the two 3D integrations over the holes momenta (hi,h 2 ) by 
fixing them to some, in principle unknown, values inside the Fermi sphere, while keeping the 
integration over the emission angle and constraining it by energy conservation. Of course, the 


obvious particularization of the mean-value theorem (2.7) is: 


F(cu,^) = Id\id^h2d^p\ /(hi,h2,p;) = Id^p\ /((hi),(h2),p;) (2.8) 


where ((hi), (h 2 )) are the two unknown hole momenta inside the Fermi sphere. 

Up to now the whole discussion is exact. What makes the difference between an excellent ap¬ 
proximation to the true result or a poor one is just the selection of the “average” hole momenta 
((hi),(h 2 )). Our choice to call it frozen nucleon approximation will be (0,0). There are mainly 
two arguments which favor this choice: 


• For high q >> kf > hi, one can assume both initial nucleons at rest. 


• If one does not assume the above statement, one has to determine the average angles for the 
two holes. And this cannot be done without computing the full 7D integral (see section 2.2 ). 

Now, we can define the frozen approximation phase-space function E{co,q)\ 

F{co,q) = {^Tte^ jd^p\5{E[+E'2-co-2m2,)&{p\,p'2,0,0)^^ (2.9) 
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where now P 2 = q — P i ■ 

The above phase-space function ( |2.9| ) allows us to define the angular distribution in the frozen 
approximation: 


F{(0,q) = ( 


2n / dd[^{d[) 

Jo 


> 


u 


> 

CD 
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d>(eO = sine; dp[p'fd{E[+E!^-(o-2mN)e{p\,P2A0) 
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( 2 . 10 ) 

( 2 . 11 ) 


Figure 1: Angular distributions 4>(0() for q = 500 MeV/c (left panel) and for ^ = 3 GeV/c (right 
panel) in the frozen nucleon approximation, for three different values of ft) as a function of the 
emission angle . 


The problem stands on the divergence of the angular distribution for some kinematics. This can be 
seen in some of the panels of figure || for fhe frozen nucleon approximafion. For ofher values of 
the holes momenta, the position of the pole is different, but it is still present. This makes crucial 
to determine analytically the angular interval where the integration has to be performed. This was 
done explicitly on section VI of Ref. []I]]. 


Once the problem of the divergence has been correctly addressed, the results for the phase 
space function in the frozen approximation can be compared with the full integral in figure 
Here we can appreciate the quality of the approximated result for a wide range of momentum 
transfers. The main discrepancies arise on the low energy-transfer region in each plot. But the phase 
space function is well-reproduced in the rest of the interval and, furthermore, we have skipped 6 
additional integrals (over holes momenta), thus reducing significantly the computation time. These 
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Figure 2: Comparison between frozen nucleon approximation and full integral for low and inter¬ 
mediate momentum transfers (left panel) and for high momentum transfers (right panel). 


results indicate that this approximation is especially well-suited for Montecarlo event generators, 
particularly if the goal is to estimate quite accurately the cross section for multi-nucleon knockout 
in the shortest time as possible. 

2.2 Other initial configurations 

We have chosen up to now the. frozen nucleon approximation, but other initial ((hi),(h 2 )) 
configurations could have been selected as well. We have considered six configurations depicted 
in fig. ^ The configurations wifh fofal inifial momenfum of fhe pair equal fo zero (U,D or T,—T) 
would, in principle, yield a resulf very close fo fhaf of the frozen approximation. This can be easily 
observed in figure ^ 

In fhe firsf comparison. Fig. we show fhe confribufion of fwo pairs of nucleons wifh fhe 
same momenfum hi = /12 = 200 MeV/c, and bofh parallel, pointing upward (U) and downward (D) 
wifh respecf fo fhe z axis, fhaf is, fhe direcfion of q. The confribufion of fhe UU configuration is 
smaller fhan average, while fhe DD is larger. This is so because in fhe UU case fhe fofal momen¬ 
fum p' in fhe final sfafe is large. By momenfum conservafion, fhe momenfa p\ and musf also be 
large. Therefore, fhese sfafes need a large excifafion energy, and fhey sfarf fo confribufe for high ft) 
fransfer. In fhe DD configurafion, fhe fotal momenfum p' is smaller, so fhe final momenfa p[ and 
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Figure 3: Geometry employed for emission of a pair of nucleons with initial parallel momenta 
(UU,DD), anti-parallel (UD,T-T) and perpendicular (TT’,-T-T’) on the left panel. In the right 
panel a comparison between average momentum approximation and full integral is shown. Initial 
momenta are both 200 MeV/c. 


p '2 can also be small, with small required excitation energy. Therefore, they start to contribute at 
lower ft). 

In the example of Fig. Ei’ two anti-parallel configurations are shown. In the UD case, one nucleon 
is moving upward and the other downward the z axis with total momentum zero of the pair. This 
situation is similar to that of a pair of highly correlated nucleons with large relative momentum 
[]_3]. Since the total momentum is zero, the final 2p-2h sfafe has fofal momenfum q, exacfly fhe 
same fhaf if would have in the. frozen nucleon approximation. Therefore, fhe confribufion of fhis 
configuration is similar to the average. The same conclusions can be drawn in the case of the con¬ 
figuration T, —T, with one nucleon moving along the x axis (transverse direction) and the other 
along —X with opposite momentum. The contribution of this pair is exactly the same as that of the 
UD configuration in the total phase-space function. 

Finally, we show in Fig. El two intermediate cases that are neither parallel nor anti-parallel con¬ 
figurations. They consist of two pairs of transverse nucleons moving along mutually perpendicular' 
directions. In the first case, we consider a T nucleon and a second T' nucleon moving in the y axis 
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Figure 4: Comparison between average momentum approximation and full integral. In the left 
panel the initial momenta are both 200 MeV/c with total momentum of the pair equal to zero. In 
the right one the initial momenta are both 200 MeV/c pointing in orthogonal directions. 


out of the scattering plane. The contribution of the TT' pair is large, while the one of the opposite 
case, —T, —T', is small. On the average, they are close to the total result. 

3. Conclusions 

We have performed a detailed study of the two-paiticle- two-hole phase-space function, which 
is proportional to the nuclear two-particle emission response function for constant current matrix 
elements. Our final goal was to obtain accurate enough results without calculating the 7D inte¬ 
gral. The frozen nucleon approximation (ID integral), that is, neglecting the momenta of the initial 
nucleons for high momentum transfer, seems to be a quite promising approach to reduce the com¬ 
putation time without missing significant accuracy. The CPU time of the 7D integral has been 
reduced significantly. We are presently working on an implementation of the present method with 
a complete model of the MFC operators. 
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